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^ ; ABSTRACT 

>J ■ We study the relation between molecular gas and star formation in a volume-limited 

sample of 222 galaxies from the COLD GASS survey, with measurements of the CO(l- 
0) line from the IRAM 30m telescope. The galaxies are at redshifts 0.025 < z < 0.05 
and have stellar masses in the range 10.0 < log M*/M@ < 11.5. The IRAM measure- 
ments are complemented by deep Arecibo HI observations and homogeneous SDSS 
and GALEX photometry. A reference sample that includes both UV and far-IR data 
is used to calibrate our estimates of star formation rates from the seven optical/UV 
] bands. The mean molecular gas depletion timescale (id ep (H2)) for all the galaxies 

in our sample is 1 Gyr, however ^^(Eb) increases by a factor of 6 from a value of 
~ 0.5 Gyr for galaxies with stellar masses of ~ 10 10 M Q to ~ 3 Gyr for galaxies with 
masses of a few xlO 11 M0. In contrast, the atomic gas depletion timescale remains 
contant at a value of around 3 Gyr. This implies that in high mass galaxies, molec- 
ular and atomic gas depletion timescales are comparable, but in low mass galaxies, 
molecular gas is being consumed much more quickly than atomic gas. The strongest 
dependences of idej^Eb) are on the stellar mass of the galaxy (parameterized as log 
£d e p(H 2 )= (0.36 ± 0.07)(logM* - 10.70) + (9.03 ± 0.99)), and on the specific star 
formation rate. A single td ep {^-2) versus sSFR relation is able to fit both "normal" 
star-forming galaxies in our COLD GASS sample, as well as more extreme starburst 
galaxies (LIRGs and ULIRGs), which have td ep {)^2) < 10 s yr. Normal galaxies at 
z—1-2 are displaced with respect to the local galaxy population in the idepCH^) versus 
sSFR plane and have molecular gas depletion times that are a factor of 3-5 times 
longer at a given value of sSFR due to their significantly larger gas fractions. 

Key words: galaxies: fundamental parameters - galaxies: evolution - galaxies: ISM 
- radio lines: galaxies - surveys 
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1 INTRODUCTION 

The physics governing the formation of stars is highly non- 
linear and operates over a vast range in physical scale. 
The size of the self-gravitating cores of molecular clouds 
in which individual stars form is only ~ 0.1 pc and the 
evolution of these cores is governed by an array of highly 
complex physical processes: turbulence and magnetic fields 
in the interstellar medium may regulate the rate of col- 
lapse of these structures, while winds, jet-like outflows, radi- 
ation pressure and ionizing radiation from the most massive 
stars will act to suppres s the formation of young stars (see 
iMcKee fc Ostrikerll2007l . for a recent review). 

It is therefore not surprising that the rate at which gas 
is being converted to stars in a single galaxy exhibits large 
scatter fro m one region to ano ther when studied on small 
scales (e.g. ISchruba et al]|2010f) . What is perhaps more re- 
markable is the fact that once the smoothing scale exceeds 
a radius of ~ 1 kpc, well-defined relations between star for- 
mation and the global gaseous properties of galaxies begin 
to appear. 

The fir st study of larg e-scale star formation in galax- 
ies was by ISchmidtl (Il959l ). who studied the relative dis- 
tribution of atomic gas and young stars perpendicular to 
the galactic plane, and derived a power-law relation be- 
tween the rate of star formation and the volume density 
of interstellar gas with exponent n ~ 2. Almost all stud- 
ies since then have compared the surface density of newly- 
formed stars with that of the total (atomic+molecular) 
cold gas (e.g lKermicuttlll989l . Il998bl : iKennicutt et al.ll2007l ; 
IBigiel et alj|2008h . These studies have confirmed the exis- 
tence of a power-law relation of the form Esfr = aE^ s , 
with N generally in the range 1.3.-1.4. 

It should be noted, however, that the exponent in the 
power law is quite sensitive to the gaseous tracer that is used, 
as well as to the type of galaxy under investigation and the 
size scale being considered. The steepest relations are gen- 
erally obtained for lower mass disk galaxies in wh ich atomic 
hydrogen dominates th e total gas surface density (|Kennicuttl 
ll989l ; lL~erov et al1l2008h . In the inner disks of the more mas- 
sive nearby spirals, molecular gas dominates. In recent work 
that made use of high-resolutio n HI, CO, 2 4 um a nd UV 
maps of nearby spiral galaxies, IBigiel et all (|2008h found 
that Esfr, and Eh2 were related by a power law relation with 
N = 1, i.e. they found that H2 appeared to form stars at a 
constant efficiency. The range of Eh2 over which this relation 
was applicable was between 3 and 50 Mq pc -1 , i.e. signifi- 
cantly lower than gas densities in starburst galaxies. In sam- 
ples of infrared-selected galaxies, including ultra-luminous 
infrared galaxies with Lir > W 12 Lq, there is an almost 
linear relation between star formation rate and the total 
amount of ver y dense molecular gas t r aced by HCN rather 
than by CO (|Gao fc Solomon! 12004 | Gracia-Carpio et al.l 
2008). Taken together, these results may imply that the 
amount of dense gas available for star formation is higher 
in more actively star-forming galaxies, but that SFR in the 
very dens est phases may be indep endent of global galaxy 
poperties (|McKee fc Ostrikeij[2007l '). 

Most of these results are, however, based on the study 
of single or of relatively small samples of objects with 
similar properties (either nearby star-forming spiral galax- 
ies or IR-selected samples of starbursting systems). Re- 



cently, ISchiminovich et al.l (|2010l ) studied the relation be- 
tween atomic gas content and star formation using an unbi- 
ased volume- limited sample of galaxies with M* > 1O 1O M0 
from the GALEX Arecibo SDSS Survey (GASS). They 
found an average HI gas consumption timescale of ~ 3 Gyr. 
Although the galaxy-to-galaxy scatter around this average 
value was large, they did not find any dependence of the ef- 
ficiency with which HI is converted into stars on parameters 
such as stellar mass, stellar mass surface density, colour, and 
concentration index. This result is perhaps not surprising, 
since the atomic gas is one step removed from the actual 
star formation process. 

The relation between atomic gas, H2 and star 
formation was also extensively investigated by the 
THINGS /HERACLES team in a series of papers 
i Walter et al l 120081 ; iLerov et al l |200S| ; IBigiel et al.fl2008l ; 
iLerov et alj l2009h . In their sample, the HI-to-I/2 ratio 
varies with radius, stellar surface density and pressure. 
These authors suggested that mechanisms occ urring at the 
GMC scale were responsible for these trends (|Lerov et al.l 
2008). Another possible scenario is that some fraction of 
the atomic gas in galaxies should be regarded as repository 
of baryonic material that has been accreted from the 
surrounding intergalactic medium. Theory predicts that 
this material re sides predomina ntly in the outer regions 
of galaxies (e.g. IMo et all Il9981 ) and much of it therefore 
has not yet reached high enough densities for self-shielding 
from the ambie nt UV field to permit molecule formation to 
take place (e.g. IFu et ai]|2010l '). 

Once molecular gas is formed, however, the picture 
painted by the THINGS/HERACLES analyses of normal 
spirals is very simple: one would expect star formation to 
occur at a rate that would deplete the H2 in a universa l 
timescale of ~ 2 Gyr (|Lerov et al.ll2008l ; IBigiel et al] 120081 ). 
Intriguingly, the same universal timescale also appears to 
hold for the "normal" pop ulation of star-form ing disk galax- 
ies at redshifts 1 < z < 2 (|Genzel et al.ll2010l ). 

We are conducting the COLD GASS survey, which will 
obtain CO (1-0) line measurements and hence the molec- 
ular gas content of a sample of ~ 350 massive galaxies 
(10.0 < log M»/iW© < 11.5) using the IRAM 30m telescope. 
The details of the survey, data for the first h alf of the sample, 
and fi rst science results were presented in ISaintonge et al.l 
Hon]) (hereafter, Paper I). Since the COLD GASS sam- 
ple is selected only by stellar mass, it enables us to study 
the relation between molecular gas and star formation for a 
representative sample of galaxies, and to test whether the 
universal depletion timescale that is observed to hold when 
the data is smoothed on scales of a few hundred parsecs in 
small samples of nearby spirals, is truly applicable to the 
galaxy population at large. 

In Section[2]we describe the COLD GASS sample of lo- 
cal massive galaxies and present the different data sets and 
measurements used in this paper. In Section[3]we present our 
main result concerning the non-universality of the molecu- 
lar gas depletion timescale, which we summarize in Section 
[S]and discuss in Section [S] Finally, in Section [3] we compare 
our results for the large and unbiased COLD GASS sample 
with results from the literature on molecular gas depletion 
timescales in normal and starburst galaxies, both at low and 
at high redshifts. Throughout the paper, distance-dependent 
quantities are calculated for a standard flat ACDM cosmol- 
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ogy with Hq — 70km s _1 Mpc -1 , and we adopt a conversion 
factor from CO luminosity to H2 mass of aco = 3.2M© (K 
km s" 1 pc 2 ) -1 (which does not account for the presence of 
Helium), unless otherwise specified. 



2 DATA 

The COLD GASS sample is stellar mass-selected (M, > 
10 10 M Q ) and volume-limited (0.025 < z < 0.050). The 
galaxies are selected from the area of sky covered by the 
SDSS, the GALEX Medium Imaging Survey (MIS) and the 
ALFALFA HI survey. This provides us with optical imaging 
and spectroscopy (over the central 3") for all galaxies in our 
sample, as well as HI fluxes for the most gas-rich systems. 
For the remaining galaxies, the HI data is obtained through 
pointed observatio ns at the Arecibo obs ervatory as part of 
the GASS survey |Catinella et al.ll2010T i. The sample used 
here contains the 222 galaxies presented in Paper I. The 
sample selection, observations and data products are fully 
described in that paper, but we present a short recapitula- 
tion here. 



2.1 Optical and UV photometry 

Parameters such as redshifts, sizes, magnitudes, and Galac- 
tic extinc tion factors are retrie ved from the SDSS DR7 
database (|Abazaiian et al]|2009T ). The UV data are taken 
from the GALEX All-sky a nd Medium Imaging surveys ( AIS 
and MIS, respectively, see iMartin et al.l 120051 ). The SDSS 
and GALEX ima ges are, however, reprocessed following 
IWang et al.l (|2010h . in order to obtain accurate UV/optical 
matched-aperture photometry. This process includes ac- 
curately registering the GALEX and SDSS images and 
smoothing them so that they have the same resolution. In 
practice, the SDSS r— band images are convolved to the reso- 
lution of the GALEX UV images before Sextractor is used to 
calculate magnitudes in consistent apertures, therefore en- 
suring that measurements in different bands cover the same 
physical regions of the galaxies. The derived NUV— r colours 
ar e corrected for Gala ctic extin ction using the presc ription 
of lWvder etafl l|2007l ) (see also ICatinella et alj|2010h . Stel- 
lar masses are ca lculated from the S DSS photometry using 
the technique of ISalim et al.l |2007l ) assuming a Chabrier 
IMF. The optical- and UV-derived quantities used through- 
out this study are published in Table 1 of Paper I. 

In this study, we also consider morphological par ame- 
ters derived for the GAS S sample by Wang et al.l (2010), fol- 



lowing the technique of lLotz et al. |2004 ). The asymmetry 



parameter (Asym) provides a measure of the difference be- 
tween the SDSS g— band image of a galaxy and a copy of the 
same image, rotated by 180°. The larger the value of Asym, 
the more asymmetric the galaxy. The other parameter we 
investigate in our analysis is smoothness (Smth), computed 
from the difference between the original g— band image and 
a copy of the same image smoothed to a scale of 20% of 
the petrosian radius. Galaxies with clumpier morphologies 
therefore have larger values of Smth. 



2.2 Optical spectroscopy 

For 113 of the COLD GASS galaxies in the present sample, 
long-slit optical spectra were obtained with the MMT 6.5m 
and the Apache Point Observatory (APO) 3.5 m telescopes. 
The spectra were taken along the major axis of the galax- 
ies and cover the wavelength range 3900-7000 A. Data re- 
duction includes standard procedures: biasing, flat-fielding, 
cosmic ray rejection, background estimation, and finally reg- 
istration and co-addition of the individual frames. Flux cal- 
ibration is achieved through observations of spectrophoto- 
metric standard stars and by bootstraping to SDSS pho- 
tometry measured through equivalent apertures. These final 
calibrated sp ectra are then fitted with a family of templates 
derived from iBruzual fc Charlotl |2003t ) stellar population 
models following t h e app roach of Tremonti et al (2004) and 
iBrinchmann et al] |2004r ). After subtraction of the stellar 
continuum, the strength of key emission features is measured 
by Gaussian fitting and the emission line fluxes are corrected 
for dust extinction using the measured Balme r decrement. 
These procedures are described in detail in IMoran et al.l 



2.3 Arecibo observations 

The HI obs e rvatio ns are described in detail in 
ICatinella etHl (|2O10Tl . so we only provide a brief overview 
here. The GASS survey builds u pon existing HI datab ases: 
the Cornell digital HI a rchive JSpringob et all 120051 ') and 

305fl 



the ALFALFA survey l|Giovanelli et all 120051 ). HI data 
for about 20% of the GASS sample (the most gas-rich 
objects), can be found in either of these sources. For the 
rest of the sample, observations are carried out at the 
Arecibo Observatory. Integration times are set so as to 
detect HI gas mass fractions (/hi —Mm/M*) of 1.5% or 
more. Observations are carried out using the L— band Wide 
receiver and the interim correlator, providing coverage of 
the full frequency interval of the GASS targets at a velocity 
resolution of 1.4 km s _1 before smoothing. Data reduction 
includes Hanning smoothing, bandpass subtraction, radio 
frequency interference (REI) excision, flux calibration and 
weighted combination of individual spectra. Total Hl-line 
fluxes, velocity widths and recessional velocities are then 
mea sured using linear fitting of the edges of th e HI profiles 
(e.g. ISpringob et al.ll2005l : ICatinella et al.ll2007l ). 



2.4 IRAM 30m observations 

We measure the CO(1-0) rotational transition line using the 
IRAM 30m telescope and its Eight Mixer Receiver (EMIR). 
As a backend, we use the Wideband Line Multiple Autocor- 
relator (WILMA), giving a resolution of 2 MHz (~ 5 km s _1 
at the observing frequency). We also simultaneously record 
the data with the 4MHz Filterbank, as a backup. With a sin- 
gle tuning of the receiver at a frequency of 111.4081 GHz, 
we are able to detect the redshifted CO (1-0) line for all the 
galaxies in our sample (0.025 < z < 0.05), within the 4 GHz 
bandwidth covered by the spectrometers. We integrate on 
each galaxy until the first of the following three conditions 
is met: (1) the CO line is detected with S/N > 5, (2) sen- 
sitivity to a molecular gas mass fraction (Mm/A/*) of 1.5% 
is achieved, or (3) an absolute minimum rms of l.lmK (per 
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20km s _1 -wide channel) is reached. For most galaxies, the 
total CO line flux can be recovered by a single pointing of 
the IRAM 30m telescope. For the most extended objects, a 
second pointing is performed, offset by 16" along the ma- 
jor axis. A set of models is used to apply the appropriate 
aperture corrections based on these measurements. 

The data are reduced with the CLASS/GILDA93 soft- 
ware. All scans are visually examined, and those with dis- 
torted baselines, increased noise due to poor atmospheric 
conditions, or anomalous features are discarded. The indi- 
vidual scans for a single galaxy are baseline-subtracted (first 
order fit) and then combined. This averaged spectrum is fi- 
nally binned to a resolution of ~ 20km s _1 . 

The flux in the CO (1-0) line is measured by adding the 
signal within an appropriately defined windowing function. 
If the line is detected, the window is set by hand to match 
the observed line profile. If the CO line is undetected or very 
weak, the window is set either to the full width of the HI 
line (W50hi) or to a width of 300 km s" 1 in case of an 
HI non-detection. The detection rate of the CO line is 54%, 
but our observing strategy allows us to set a stringent upper 
limit in the case of a non detection. 

After correcting the CO (1-0) line fluxes (or upper- 
limits) for aperture effects, using a set of models de- 
rived from a compilation of nearby galaxies with resolved 
maps, total CO lum inosities L' co are calculated following 
ISolomon et al.l (|l997T ). The total molecular hydrogen masses 
are then calculated as Mh2 = L' co aco- We adopt a con- 
stant Galactic conversion factor of aco = 3.2 Mq (K km s _1 
pc 2 ) -1 , which does not include a correction for the presence 
of Helium. 

An extensive description of the IRAM observing proce- 
dure and of the data reduction steps can be found in Paper 
I. The CO line fluxes and molecular gas masses for all the 
galaxies in this study are given in Table 2 of Paper I. 




8 9 10 11 

log L co [K km s -1 pc 2 ] 

Figure 1. Relation between CO luminosity and far-infrared 
(FIR) luminosity. Small gray circles show the compila- 
tion of normal and lumi nous infrared galaxies presented in 
iGracia-Carpio et al. The dashed line shows the relation 

derived by Ge nzel et al for normal galaxies (the shaded 

region denoted the Icr scatter about this relation) ;the dotted line 
is the relation derived by the same authors for ULIRGs/mcrging 
systems. COLD CASS galaxies with IRAS measurements are 
overplotted as filled circles and are seen to follow the same trend. 
For the remainder of the COLD CASS sample (open blue cir- 
cles), we "infer" Lfir from the total star formation rate after 
subtracting the component contributed by the observed UV ra- 
diation (without attenuation correction). 



2.5 Star formation rates 

Star formation rates for all COLD CASS galaxies are calcu- 
lated from optical/UV spectral energy distribution (SED) 
fitting (SFRsbd). The technique is explained in detail in 
Appendix[A] In short, fluxes in the seven SDSS and GALEX 
bands are fitted with a serie s of mo del SEDs constructed 
from the iBruzual fc Chariot! ||2003T ) population synthesis 
code, assuming a range of star formation histories, ages, 
metallicities and dust attenuation factors. As described in 
the Appendix, the assumed distribution of dust attenua- 
tion factors is the critical factor in obtaining accurate star 
formation rates usin g this method. We use the sample of 
Ijohnson et al.l {2007]) to calibrate our dust attenuation pri- 
ors. This dataset includes GALEX FUV data, plus 8, 24 
and 70 /im data from Spitzer for a representative sample of 
galaxies selected from SDSS. For galaxies with NUV— r< 5, 
our values of SFRsed are consistent with the total star for- 
mation rates (SFRtot) derived from the combined GALEX 
FUV and Spitzer 70/im fluxes, with scatter of only 0.22 dex 
(see Appendix [Al and Figure IAT1) . 



1 |http : //www . iram.fr /IRAMFR/GILDAS 



2.6 Comparison with infrared luminosity 

For a small fraction of the COLD GASS sample (~ 15%), 
flux measurements at 60 and 1 00 fim can be fou nd in 
the IRAS Faint Source Catalog (jMoshir et all 1 19921 ). The 
far-infrared luminosity inferred from these fluxes is an ex- 
cellent measurement of the total luminosity from young 
stars that has been reprocessed by dust, and is known 
to correlate well with the total CO luminosity of the 
galaxies (e.g. Sanders fc Mirabel 1 19851 : ISanders et al.1 Il99ll ; 
iGao fc Solomonll2004f ). In Figure IT1 we plot the Lfir — Lco 
relation for the COLD GASS sample and compare it to liter- 
ature results for a collection of nearby galaxies. The COLD 
GASS galaxies with IRAS measurements (filled circles in 
Figure [1} follow previously derived relations almost exactly. 
The rms scatter about the relation is 0.30 dex, compared to 
the value of 0.31 dex derived for reference samples of nearby 
galaxies l|Genzel et alj|20ich . 

For the COLD GASS galaxies without IRAS measure- 
ments, we infer Lfir by assuming that it should be equal 
to (SFRsED-SFRfyv")- These estimates are plotted as open 
sy mbols in Figu r e CD F or the sample, the scatter around 
the lGenzel et al.1 (|2010h relation is 0.39 dex. This enhanced 
scatter could be due to the uncertainty in the values of 
SFRs_B_o(see Appendix [XJ; only direct infrared measure- 
ments of these galaxies will allow us to firmly establish the 
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Figure 2. Molecular gas mass fraction as a function of NUV-r 
colour for the current COLD GASS sample of 222 galaxies. Secure 
detections are represented as filled circles, tentative detections as 
open circles, and non-detections as 5<r upper limits as indicated 
by the arrows. The best fit linear regression is overplotted, as well 
as the empirical threshold of NUV— r= 5 which marks an abrupt 
transition from the gas-rich, star forming population to a more 
passive population. 



2.8 Sample definition 

One of the key results of Paper I is presented in Figure [2] We 
found sharp thresholds in galaxy parameters such as stellar 
mass (M t ), stellar mass surface density (/x*), concentration 
index (R90/R50) and colour, below which almost all galaxies 
have a measurable cold gas component but above which the 
detection rate of the CO line drops suddenly. These thresh- 
olds correspond approximately to the transition between 
blue cloud and red sequence. The strongest bimodality is 
seen when relating the molecular gas mass fraction and the 
NUV— r colour of the galaxies (Figure [2]). A clear threshold 
at NUV— r~5 marks a break between the star forming and 
the passive populations. 

This colour threshold is also the limit where the un- 
certainty on the SED fitting-derived star formation rates 
increases abruptly (Figure [Al[) . Galaxies with NUV— r> 5 
only have upper limits for Mh 2 and very uncertain measures 
of SFR, therefore their depletion times (Equation]!]) are very 
poorly constrained. From here on, we therefore consider only 
the subset of COLD GASS galaxies which have NUV-r< 5 
and which make up the active population. Through stack- 
ing, we have shown in Paper I that the passive population 
is truly "red and dead" , with a stringent constraint on their 
average gas mass fraction of Mh 2 /M, — 0.0016 ± 0.0005. 
These red galaxies are therefore mostly irrelevant for the 
following discussion and we do not consider them further. 
Galaxies with NUV— r< 5 but no CO detection are consid- 
ered as upper limits in the tdepiH?) plots, using 5a upper 
limits for Mh 2 (see Paper I). 



level of scatter in the Leo — Lfir relation for our com- 
plete sample. Nevertheless, Figure [1] makes it clear that our 
sample spans a broader range in properties such as star for- 
mation rate compared to IRAS-selected samples. As we will 
show in Section [3l our samples also include galaxies with a 
wider range of CO-to-IR luminosities or, equivalently, molec- 
ular gas mass to star formation rate ratios. 



3 A NON-UNIVERSAL MOLECULAR GAS 
DEPLETION TIMESCALE 

Resolved studies of the disks of normal galaxies indicate 
that the molecular gas depletion timescale is constant 
l|Bigiel et a l. 2008; Lerov et alj|2008h . In this section, we in- 
vestigate how our integrated measurements of tdep(H2) scale 
with a number of other global physical parameters, including 
stellar mass, stellar surface mass surface density, bulge-to- 
disk ratio and specific star formation rate. 



2.7 Molecular gas depletion timescales 

Using the molecular hydrogen gas masses (or upper limits) 
and the star formation rates derived from SED fitting, we 
define for each galaxy its molecular gas depletion timescale: 

t dep (U 2 ) = Mh 2 SFRsed~\ (1) 

This parameter therefore describes how long each galaxy 
could sustain star formation at the current rate before 
running out of fuel, assuming that the gas reservoir is 
not replenished. Some studies instead consider the inverse 
of tdep (H2), the molecular gas star formation efficiency 
(SFE_h 2 ). However, since the definition of SFE is not uni- 
versal, with some authors quoting SFE as a dimension- 
less quantity per dynamical time, we use in this study the 
more straightforwardly-defined depletion timescale, defined 
in Equation [T] Similarly, the atomic gas depletion timescale 
is t dep (m)= MhiSFRsed" 1 . 



3.1 Variations with global physical parameters 

In Figure [3] the molecular gas depletion timescale is plotted 
against a range of parameters describing the global physical 
properties of the galaxies. Although there is scatter, a clear 
dependence of ide^H^) on many of the parameters is ob- 
served. In all cases, we use two different measures to quantify 
the strength of the dependence of id ep (H2) on the x-axis pa- 
rameter: (1) the Pearson correlation coefficient of the param- 
eters, r, and (2) the scatter of the residuals in log(td ep (rl2)) 
of the best fit linear relation (regression of y against x), a. As 
a comparison point, the scatter in log(f<j e p(H2)) for all the 
CO detections is a = 0.315 dex. We fit both to the secure de- 
tections only, and also to the entire sample of galaxies with 
NUV— r< 5 including tentative detections and upper lim- 
its. The coefficients as well as the parameters of the best-fit 
relations for both of these (sub-)samples are given in Table 
[1] All the fits and mean values are computed by weighting 
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Figure 3. Depletion time for the molecular gas as a function of a variety of parameters. Filled and open black circles indicate 
COLD GASS secure and tentative detections. We also show the mean values in equally populated bins for secure detections 
only (blue symbols) and for the entire sample with NUV— r< 5 (red symbols). The error bars indicate la confidence interval 
on the mean, determined by bootstraping. The solid and dashed lines are the linear regression of ttjep(H2) against the x— axis 
parameter, for the secure detections only and the entire sample, respectively. In each panel we give the scatter in the residuals 
of this best- fit relation (a), as well as the correlation coefficient measured for the secure detections (r). 



galaxies differently according to their stellar mass. This pro- 
cedure, explained in detail in Paper I, corrects for the fact 
that by selection the COLD GASS sample has a log M* dis- 
tribution which is flatter than in a purely volume-limited 
sample. 

Interestingly, the depletion timescale depends both on 
parameters that describe the masses and structural prop- 
erties of the galaxies and on parameters relating to their 
stellar populations and star formation rates. A clear de- 
pendence of tde P (H2) on the stellar mass of the systems 
is observed (r = 0.48), as well as on the presence of a 
bulge as parametrized by the concentration index, -R90/-R50 
(r = 0.40), and on pressure within the disk as measured by 



the stellar mass surface density, (r = 0.44). In all these 
cases, the scatter in the residuals of the best-fit relations is 
significantly smaller than the scatter in tdepdHb) itself. Of 
these three, the dependence on log M» is strongest, and also 
the most insensitive to the presence of non-detections in CO 
(see Figure^ and Table (TJ. 

There is a weaker dependence of td ep (R2) on dynamical 
mass (Figure [3jl, r = 0.31), which is estimated as: 

M dyn = (2) 

where the circular velocities, v c , are measured from the HI 
linewidths, W50hi, as: 
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Table 1. Molecular gas depletion time relations' 2 



All galaxies with NUV — r < 5 Secure detections only 



x parameter units xq m bar m bar 



log M* 


log M e 


10.70 


0.36 ± 0.07 


9.03 ± 0.99 


0.274 


0.42 


0.40 ± 0.07 


9.06 ±1.03 


0.269 


0.48 


log (i. 


log M Q kpc -2 


8.70 


0.33 ± 0.07 


8.95 ± 0.91 


0.284 


0.36 


0.40 ± 0.07 


8.99 ± 0.88 


0.275 


0.44 


log SFR/M, 


log yr -1 


-10.40 


-0.44 ± 0.04 


8.98 ± 0.61 


0.246 


-0.61 


-0.54 ± 0.04 


9.06 ± 0.62 


0.221 


-0.70 


log M dyn 


logMg 


11.40 


0.28 ± 0.07 


8.98 ± 1.05 


0.293 


0.31 


0.27 ± 0.08 


8.99 ± 1.16 


0.291 


0.31 




2.50 


0.30 ± 0.06 


8.92 ± 0.21 


0.286 


0.34 


0.37 ± 0.07 


8.99 ± 0.22 


0.290 


0.40 


log E SFR 


log M yr _1 kpc _2 


-2.40 


-0.24 ± 0.04 


8.97 ± 0.15 


0.283 


-0.43 


-0.24 ± 0.06 


9.01 ± 0.16 


0.287 


-0.44 


NUV-r 


mag 


3.50 


0.16 ± 0.03 


8.94 ± 0.14 


0.275 


0.42 


0.18 ± 0.03 


8.98 ± 0.15 


0.263 


0.49 


Smth 




0.05 


0.46 ± 1.01 


8.93 ± 0.08 


0.304 


0.06 


1.60 ± 1.28 


8.93 ± 0.10 


0.307 


0.13 


Asym 




0.10 


0.05 ± 0.49 


8.93 ± 0.08 


0.305 


0.10 


0.34 ± 0.61 


8.94 ± 0.09 


0.309 


0.16 



a The relations are parametrized as log t<jep(H2)[yr x ] = m(x — xq) + b, with all quantities and units given in this table. 



WS0hi/2.0 

Vc - (l + z)sin(incl)' [ ) 

The factor (1 + z) accounts for cosmological stretching and 
the sin(incl) term is a correction for the inclination of the 
system (as measured from SDSS photometry). In Equation 
[2] we further assume that the size of the HI disk (7?) is 1.5 
times the size of the optical disk (D25/2). Uncertainty in 
this approximation certainly contributes to the scatter of 
the relation between td ep (R2) and Md yn - 

The strongest dependencies in Figure [3] are however on 
quantities that are sensitive to the amount of star formation 
in the galaxies. The two strongest correlations of td ep {^2) 
are with NUV— r colour (r = 0.49) and with specific star 
formation rate (sSFR, r — —0.7). These quantities are di- 
rectly linked, because NUV— r is a good proxy for sSFR since 
it relates a quantity tracing ongoing star formation activity 
(NUV flux) and a quantity sensitive to the older stellar pop- 
ulation (r— band flux). However, they differ in the sense that 
the sSFR takes into account internal dust attenuation and 
the mass-to-light ratio of the old stellar population, while 
NUV— r is not corrected for both these effects. 

At first glance, id e p (H2) seems to depend more strongly 
on sSFR than on NUV— r colour, as revealed by the corre- 
lation coefficients and scatter around the best fit relations 
(Table [T]). However, we must be cautious because the same 
measure of the SFR enters in the computation of tde P (R2) 
and sSFR, and the stronger correlation with sSFR compared 
to NUV— r may therefore only be induced by this common 
quantity. We test for this using a Monte-Carlo approach, 
which aims to decorrelate the errors on td ep (H-2) and sSFR. 
To mimic an independent measure of the two parameters, 
while keeping tde P (H-2) fixed to the original values, we re- 
compute the sSFR by allowing the star formation rate of 
each galaxy to vary around the measured value according to 
the observed measurement error, assuming a normal distri- 
bution (asFR = 0.24 dex, see Figure Ell . We perform 5000 
iterations and examine the distribution of the correlation 
coefficients in the simulated relations, which has a median 
value of r = —0.61 as shown in Figure [3] While there is 
only a 0.5% chance that the original strength of the corre- 
lation (r = —0.71) is not due to measurement correlation, 
the probability that the relation is more strongly correlated 
than the i dep (H 2 )-NUV-r (r = 0.49) or the t d ep(H 2 )-M„ 
(r = 0.48) relations is > 99.9%. We therefore conclude that 
while part of the strength of the correlation between i dep (H2) 
and sSFR is indeed induced by the common measurement of 
SFR entering both parameters, the relation is still at least if 



not more strongly correlated than the next most significant 
correlations (i.e. those on NUV— r colour and stellar mass). 
In Section T3.2I below, we present further evidence for this. 

The molecular gas depletion time is also correlated 
with the star formation surface density, measured within 
-Rso(SFR), the radius enclosing half of the total star forma- 
tion measured for the whole galaxy. This quantity is mea- 
sured as part of the SFRsbd calculation pipeline. The gas 
surface density, E sas , cannot be accurately calculated for the 
COLD GASS sample since no direct measurement of the size 
of the gas disks is available. However, since E ffas and Esfr 
are directly linked through the Kennicutt- Schmidt relation, 
we can infer from Figure (3f that tde P (H-2) is a decreasing 
function of E gas . 

On the other hand, t dep (H2) does not appear to depend 
on measures of the degree of asymmetry and dumpiness 
in the galaxies (Figure [3ji,i), the scatter of the residuals of 
these two relations (a =0.308 and 0.306, respectively) is not 
significantly reduced compared to the scatter of 0.315 dex 
in tde P (H2) for the whole sample. Such dependencies might 
have been expected if, for example, spiral density waves 
or other perturbations in the disk trigger the formation of 
GMCs and thus contribute locally to a decrease in tde P (H2)- 
The fact that we do not see a clear dependene of £ dep (H2) 
on our measurements of structure within the disks could 
mean either that 5— band light (the band used to measure 
Asym and Smth) is not an optimal tracer of these compres- 
sion events, or else that the effect mostly disappears once 
averaged over entire galaxies. We note however that this 
independen ce of t deP (H2) on s piral arm patterns was also 
reported bv lLerov et alj (|2008T l using resolved data. 

To illustrate some of the results presented in Figure [3] 
graphically, we show in Figure [5] a gallery of SDSS thumb- 
nails of COLD GASS galaxies at their appropriate locations 
in the sSFR-f de p(H2) plane. The images clearly reveal that 
galaxies with longer molecular gas depletion times tend to 
be redder, more massive, and have larger bulges than those 
with shorter depletion times. 

As a final step, we investigate whether third parameter 
dependences can be found. We analyze whether the resid- 
uals of the sSFR-f dep (H2) relation (the tightest correlation 
found in Figure [3]) show any dependence on any of the pa- 
rameters included in this study (see Figure IB 1 1 and details 
in Appendix [Bj. We do not identify any convincing third 
parameter dependence. 

Our results imply that for a given amount of molecular 
gas, all galaxies are not equally efficient at processing this 
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gas into stars. We find an increase in the mean molecular gas 
depletion timescale by a factor of ~ 6 over the stellar mass 
range of 10 10 to 10 11,5 M©. This implies that galaxies with 
large stellar masses, red colours, high stellar surface densi- 
ties and/or older stellar populations will require significantly 
longer times to turn their molecular gas into stars. The mean 
depletion timescale in these systems is tde P (H2)~ 3 Gyr, as 
compared to the average of ~ 1 Gyr for the sample as a 
whole, and to the mean value of ~ 0.4 Gyr for galaxies with 
stellar masses ~ 10 10 M Q . 



3.2 Validation using spectroscopic parameters 

In this section, we analyze the behavior of tdepi^) as a 
function of two spectroscopically-derived quantities, which 
we measure either from SDSS fiber spectroscopy or from 
the long slit data (see Section f2 . 2 f) . While SDSS fiber spec- 
tra are available for the entire COLD GASS sample, they 
have the disadvantage of sampling only the central 3"of the 
galaxies. This limitation could bias our results when inves- 
tigating quantities that vary radially within the disks. 

We measu re metallicity from stro ng line ratios using the 
prescription of iPettini fc Pagell (|2004T ): 



12 + log(Q/H) = 8.73 - 0.32 log 



[Qin]/H^ 

[NH]/H Q 



(4) 



On this scale, a value of 12 + log(0/H) = 8.66 corresponds to 
solar metallicity. In Figure [BJi, we plot id ep (H2) as a function 
of metallicity measured from both sets of spectra. td ep (H2) 
is found to be completely independent of metallicity (the 
correlation coefficient is r = 0.01). In fact, of all param- 
eters investigated (e.g. in Figure [3]), metallicity is the one 
that shows the least degree of correlation with the depletion 
timescale. This result is reassuring, as a metallicity depen- 
dence could have cast some doubt on our assumption of a 
constant CO-to-H2 conversion factor (Xco)- We remind the 
reader that the COLD GASS sample only includes galaxies 
with stellar masses greater than 10 M© and less than a few 
xl0 n M©. The relation between metallicity and stellar mass 
is rather flat over this regime (Tremonti et al 2004), with the 
average gas-phase metallicity only increasing by ~ 0.25 dex 
over 1.5 dex in logM*. 

The spectroscopic data also provide us with information 
about stellar absorption features, including the strength of 
the 4000A break. The strength of this feature varies monto- 
tonically with the age of a stellar population of fixed metal- 
licity. If we assume that the galaxies in our sample have 
had smooth, exponentially declining star formation histo- 
ries, we can use standard stellar population synthesis models 
(Bruzual & Chariot 2003) to calibrate a relation between the 
4000 A break strength and sSFR. We have done this both 
for the fiber and the long-slit measurements of this feature. 

This independent estimate of the sSFR can be used as 
a second validation of the strength of the dependence of 
idep(H2) on the specific star formation rate. In Figure [Bp, 
we show that the relation between idep (H2) and sSFR es- 
timated from the 4000 A break has a correlation coeffcient 
of -0.61, almost as strong as the value of -0.70 measured for 
the same relation using photometric quantities, and iden- 
tical to the value estimated using the Monte-Carlo simula- 
tions (see Figure[4p). This spectroscopic measurement there- 
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Figure 7. Depletion time for the molecular g function 
of specific star formation rate for COLD GASS secure detec- 
tions (filled c ircles ) as well as for the HERACLES sample from 
iLerov et alj ||2008|) . We plot both the integrated data (filled 
squares) and the resolved data (crosses). The COLD GASS sam- 
ple spans a larger range in global values of sSFR than the HERA- 
CLES sample. When considering integrated measurements, both 
data sets however show similar scatter and behaviour in the range 
of sSFR where they overlap. 



fore provides independent confirmation of the effects already 
seen with NUV— r colour (Figure [3J5) and photometrically- 
derived sSFR (Figure[3p). 



3.3 Comparison with resolved studies 

We compare our results with those obtained from resolved 
studies of nearby galaxies. The two approaches are highly 
complementary: COLD GASS gives us a global picture of 
gas and star formation in the local Universe, but lacks the 
power to trace the exact distribution of these components, 
something surveys like THINGS/HERACLES are designed 
to do. The methods are however fundamentally distinct, 
since a volume-limited survey like COLD GASS gives more 
statistical weight to the lower mass galaxies which naturally 
dominate the galaxy number density, while a resolved study 
tends to give more weight to bigger, more massive galax- 
ies which cover more resolution elements. These distinctions 
should be kept in mind when comparing observations. 

While our results show a variation in tde P (H-2), in con- 
trast with the constant H2 depletion time found for HERA - 
CLES survey galaxies l|Lerov et al.ll2008l : iBigiel et al]l20Qg| ). 
the results are not surprising. HERACLES covers a limited 
range of surface densities typical of nearby "normal" spi- 
rals. Observations that probe a wider range of galaxy types 
have long been known to give different H2 depletion times 
between the regimes of normal an d "bursty" star forma - 
tion (e.g. in mergers and ULIRGS. lGao fc Solomon] l2004n . 
and theoretical models have been able to explain this di- 
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Figure 4. Validation of the t ( j ep (H2)-sSFR relation using a Monte-Carlo approach. The correlation coefficient (r) of the 
original relation is compared to the value measured for a similar relation after decorrelation of the measurement errors, (a) 
An example showing one of the repetitions and comparing the original relation (filled circles) with the "decorrelated" relation 
(open triangles), (b) Distribution of the correlation coefficient of the 5000 iterations (median r = —0.615). While there is only 
a 0.5% chance that the original strength of the correlation (r = —0.71) is not due to measurement correlation, the probability 
that the relation is more strongly correlated than the < £ j f , p (H2)-NUV— r or the t t jep(rl2)-.M» relations is 99.9%. 



chotomy as arising from a transition between galaxies where 
molecular clouds a re gravity-confined to o nes where they are 
pressure- confined <|Krumholz et al.ll2009T ). In order to illus- 
trate the consistency of our results with HERACLES, in 
Figure [7] we directly compare the two data sets. To do this 
comparison, we plot the depletion timescale as a function of 
specific star formation rate (sSFR), since this quantity can 
be reliably compared between the two samples. The compar- 
ison shows that the COLD GASS data are perfectly consis- 
tent with the HERACLES measurements, once all quanti- 
ties are brought to a common scale (same conversion factor 
a C o = 3.2M (K km s" 1 pc 2 ) -1 ). The reasons why COLD 
GASS reveals a variable molecular gas depletion timescale 
are (1) a significantly larger sample, and (2) much larger 
dynamic ra n ge in quantities such as sSFR. We note that 
iLerov et al.1 (|2008l ) quote a mean molecular gas depletion 
time of 1.9 ±0.9 Gyr. However, after using the same conver- 
sion factor and restricting their sample to the COLD GASS 
stellar mass range, the mean molecular gas depletion time 
for the HERACLES sample is 0.9 Gyr, perfectly consistent 
with the COLD GASS estimate (1.05 ± 0.74 Gyr), as seen 
in Figure [7] 

We can gain some further insight by looking at the re- 
solved data from HERACLES, also plotted in Figure [7] Not 
surprisingly, there is larger scatter in the tdep (H2)-sSFR rela- 
tion, since these points are averages over only ~kpc-sized re- 
gions and not entire galaxies. The resolved data do cover the 
full range of tdep(H2) and sSFR values probed by the COLD 
GASS galaxies, but the vast majority of the data points are, 
however, concentrated in the range (-11. 0< logSFR/M^ < 
-10.0) and (8.7 < log t dep (R 2 ) < 9.3). There are very few 
regions within these normal star forming spirals that reach 
values of td ep (H2) as low or as high as those seen in the inte- 
grated COLD GASS measurements. This suggests that large 



variations of tdep (H2) between different kpc-s ized regions 
of th e same galaxies are not common (see also ILerov et al.l 
2008), and therefore that extreme values on small scales are 
almost exclusively found in galaxies that are also globally 
extreme. The few resolved points with depletion times that 
are unusually high (low) are located in the disk centers (out- 
skirts), where the Mh2/Mhi ratio is high (low). 



3.4 Comparison with atomic gas depletion time 

Another important step in understanding cold gas deple- 
tion times in galaxies is the atomic-to-molecular transi- 
tion. In this section, we compare the molecular and atomic 
gas depletion time s cales across the COLD GASS sample. 
ISchiminovich et all (2010) studied the behavior of td ep (HI) 
for galaxies in the GASS survey. They found that the aver- 
age atomic gas depletion timescale for galaxies with M, > 
10 10 Mo is 3 Gyr. There is considerable scatter from one 
galaxy to another, but the average timescale is a constant 
function of all quantities investigated. 

In Figure [8] we compare the dependencies of tdep (H2) 
and tdep (HI) on M*, /x* and NUV— r colour. N o t su rpris- 
ingly, we reproduce well the ISchiminovich et all (|2010h re- 
sult; tdep (HI) does not depend on any of these three quanti- 
ties and is on average 3 Gyr, although the values of td ep (HI) 
for individual galaxies range from a few hundred million 
years to more than a Hubble time. In comparison, the range 
of values for tdep (H2) is smaller, from a few hundred mil- 
lion years to ~ 5 Gyr. This is not surprising, given that 
the atomic gas-to-stellar mass ratio has much larger dy- 
namic range than th e molecular gas-to-stellar mass ratio 
l|Catinella et alj|201(j Paper I). 

Using the exact same sample, we can see clearly in Fig- 
ure[8]that whereas tdep(HI) is constant at ~ 3 Gyr, tdep(H2) 
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Figure 5. SDSS three-colour images (of dimension 50x50 kpc) of a subset of the COLD GASS galaxies in the sSFR-t ( ; e p(H2) 
plane. Some of the trends present in Figure \3\ are visible: galaxies with the longest depletion times tend to be redder, with 
early- type morphologies and more prominent bulges. 



is a function of all three quantities (M», fj,* and NUV— r). 
It is interesting that at the highest masses, largest surface 
densities and reddest colours, tde P (H2)—tde P (HT) for the ma- 
jority of galaxies. However, for low mass, low stellar surface 
density and high sSFR galaxies, tde P (R2) is smaller than 
idep(HI) by almost an order of magnitude. In other words, 
low mass galaxies are much more imminently in danger of 
running out of molecular gas than out of atomic gas, whereas 
in high mass systems there appears to be a balance between 
the consumption rate of atomic and molecular gas (at least 
when averaged over entire galaxies). 



4 BRIDGING LOW AND HIGH REDSHIFT 
MOLECULAR GAS STUDIES 

While direct HI measurements in individual objects are 

still only possible for nearby galaxi es (up to z ~ 0.3, 

IVerheiien etafl 120071 ; | Catinella et al.l |2008) , the detection 
of CO line emission in galaxies at z > 1 is now achievable 
and will soon become common practice with for example 
the Atacama Large Millimeter Array (ALMA). Because of 
the unbiased and homogeneous nature of the sample, COLD 
GASS will provide an ideal local comparison point for these 
high redshift studies. In this section, we take a first look at 



how the molecular gas depletion times for the galaxies in our 
sample compare with published results at higher redshifts. 

The IRAM Plateau de Bure interferometer is up to 
now the instrument of choice to perform high redshift CO 
observations. It has been used to characterize the molec- 
ular gas content of submillimeter galaxies (SMGs, e.g 
Prayer et al] 1 19981; [5 owncs fc Solomonl 120031; iGenzel et al 



20031 ; iNeri et al.1 120031; iKneib et all 120051; [Tacconi et al 
bood : iBothwell et al.l l2010l : lEngel et al.l 120101) and norma 
star forming galax ies at z > 1 (e.g . Tacconi et al.l |2010| ; 



IGenzel et al]|2010l ; iDaddi et al.ll2010t ). In Figure [9] we com- 
pare the position of such systems in the td ep (H2)-sSFR 
plane. To perform this comparison, all data sets need to 
be put on a common scale. A standard conversion factor 
of aco = 3.2Mq (K km s _1 pc 2 ) -1 is adopted for galax- 
ies with Lfir < 10 12 Lq, and a factor of aco = I.OMq 
(K km s _1 pc 2 )" 1 is used for the ULIR Gs and merging sys- 
tems (these choices are motivated in e.g. lTacconi et ap2010l ; 
IGenzel et al.l |2010| ) . The star formation rates are derived 
from the infrared luminosity for the LIRGs and U LIRGs, 
and from the 1.4 GHz flux (|Chapman et al.l 120051 ) for the 
SMGs, using the cali bration of the FIR-radio correlation of 
iMagneUi et all (|2010l ). 

For both the high and low redshift populations, merg- 
ing systems (local ULIRGs and high-z SMGs), have much 
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Figure 6. Molecular gas depletion time as a function of two spectroscopic quantities: (a) the gas metallicity calculated using 
Equation l4l and (b) the sSFR inferred from the measurement of the 4000A break (see text). The spectroscopic quantities are 
derived from the SDSS fiber spectroscopy (small gray symbols) and independently from the long-slit spectra (large black symbols, 
open for tentative CO detections). The solid lines are the best-fitting relations to the latter. In panel (b), as a comparison, 
the la scatter around the relation of Figure [3}; derived using sSFRg^^ is shown as the shaded region. The spectroscopic data 
provide an independent confirmation of the trend observed between sSFR and t,fei,(H2). 



.5 ■ 



9.0 ■ 



60 

■2 8 



(a) 



tit-"? •* # °* 

° m e r=0.47 



9.5 



9.0 ■ 



8.5 ■ 




10.0 10.5 11.0 11.5 



8.0 8.5 9.0 9.5 




(e) 



JS>o» 



• •.. . 

• r*=-y.o 



10.5 
i 10.0 

9.5 

9.0 

8.5 



(0 



10.0 10.5 11.0 11.5 
log M. [M ] 



8.0 8.5 9.0 9.5 
log /j.. [M kpc" 2 ] 



*V * • • 

• ••• y 



• •I 



r=0.05 



3 4 
NUV-r 



Figure 8. Depletion time for the molecular gas (top row) and atomic gas (bottom row), as a function of (a,d) stellar mass, 
(b,e) stellar mass surface density and (c,f) NUV-r colour. 
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Figure 9. Comparison of the depletion time-specific star formation rate relation, at z =0, 1 and 2, for normal star forming 
galaxies and extreme systems (local LIRGs/ULIRGs and high- z SMGs). The best-fit bisector linear relation to the COLD CASS 
sample, log t dep (U2)= (—0.724 ± 0.039)sSFR + (f .54 ± 0.41), is plotted as a solid black line. This local relation is evolved to 
2=1 and 2 (d otted blue and dashed red lines, respectively) using the observed evolution in sSFR with cosmic time (using 
the relation of Perez-Gonzalez ct al. 2008). The known offset between the depletion times of normal star-forming galaxies and 
merging systems (local ULIRGs and high z submillimeter galaxies) is seen, but the population of local LIRGs is found to bridge 
this gap. A standard conversion factor of aco = 3.2Mq (K km s~ x pc 2 ) -1 is adopted for galaxies with Lpm < 10 12 Lq, 
and a factor of ago = l.OMp, (K km s —1 pc 2 ) — 1 is used for the ULIRGs and merging systems (these choices are motivated 
in e.g. iTacconi et aUl201Cl : iGenzel et alj|2010l ). The impact of using the smaller conversion factor also for the brightest LIRGs 
(Lpjjj > 10 1 1 75 L (7)) is shown by a one-side d error bar for each individual object. The high redshift data points are from 
IGenzel et ail i2010l l and lHainline et alj ll2010h , and the local LIRGs/ULIRGs are from I Howell et al.l J201(l) and lda Cunha et al.l 



shorter depletion timescales than "normal" star-forming 
galaxies. These two populations of normal and extreme 
starformers are observed as two distinct bran ches in the 
Kennicutt- Schmidt relation (IGenzel et al.ll2010l) . Usin g the 
measurements of stellar masses of IHowell etaU (|20ld ) and 
Ida Cunha et all finch for a large sample of local LIRGs, we 
see that the gap in the relation that extends from the COLD 
GASS galaxies to the ULIRGs is filled by these galaxies of 
intermediate infrared luminosities. 

We first focus our attention on the "normal" star form- 
ing population (filled symbols in Figure [9]). The z = 1 — 2 
galaxies lie at systematically larger values of sSFR as com- 
pared to our local comparison samples - this simply re- 
flects the well-studied cosmic evolution of the star forma- 



tion versus stellar mass r e lation (e.g. |Noeske et all 120071 : 
IPerez-Gonzalez et alJliooj ; iRodighiero et alj|20ld) . In Fig- 
ure [9l the f c jep(H2)-sSFR relation derived for the COLD 
GASS sample is plotted as a solid line. The dotted and 
dashed lines show this rel ation displaced along t he x- axis 
by factors derived based on lPerez-Gonzalez et al.l |2008l ) for 
the evolution of the sSFR in blue sequence galaxies out to 
redshifts of 1 and 2, respectively. Remarkably, the tde P (H2)- 
sSFR data points for the z = 1 and z — 2 galaxies lie almost 
exactly on these "shifted" relations. The reason why normal 
high redshift galaxies can have relatively long (~ 1 Gyr) 
molecular gas consumption timescales even though their sS- 
FRs have increased by nearly an order of magnitude, is be- 
cause they are so much more gas-rich than nearby objects 
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l|Tacconi et alj l2010h . The only galaxies in our local sam- 
ples that reach sSFR comparable to the normal high-z star- 
forming populations are the ULIRGs and some of the more 
extreme LIRGs, which are almost all dynamically disturbed 
or interacting systems, with much shorter molecular gas de- 
pletion times. 

Turning now to the full population of galaxies illus- 
trated in Figure [9] we see that when the COLD GASS sam- 
ple is combined with samples of low redshift LIRGs and 
ULIRGs, the idep(H2) versus sSFR relation appears to ex- 
tend (with much the same slope and even scatter) over an- 
other full 1.5 dex in log sSFR. Instead of a factor 6 range in 
molecular gas depletion time, there is a factor ~ 50 decrease 
from the quiescent end of the galaxy population where the 
present-day star formation rate is only a tenth of its past 
averaged one, to the most extreme starbursts where the 
timescale for making the entire stellar mass of the galaxy 
is 1 Gyr or less. There is a hint that the equivalent extreme 
galaxies at high-z, the SMGs, may track a different tde P (H-2) 
versus sSFR relation - one that is significantly displaced 
with respect to the relation traced by the normal starform- 
ing galaxies. Or alternatively, the slope of the tde P (H2)-sSFR 
relation may be evolving with redshift. However, testing 
whether this is the case will require further observations of 
high redshift galaxies, spanning a wider range of properties, 
especially extending to lower (specific) star formation rates. 



5 SUMMARY OF OBSERVATIONAL RESULTS 

A summary of the main empirical results presented in this 
paper is as follows: 

(i) The mean molecular gas depletion timescale for COLD 
GASS galaxies with M» > 1O 1O M and detected CO (1-0) 
line is 1.05 ±0.74 Gyr. 

(ii) We observe an increase in tde P (H2) by a factor of 6 
over the stellar mass range of 10 10 to lO n ' 5 M0, from ~ 0.5 
to ~ 3 Gyr. 

(iii) Over the same stellar mass range, the atomic gas 
depletion timescale remains contant at a value of around 3 
Gyr. This means that in high mass galaxies, the molecular 
and atomic gas depletion timescales are comparable. In low 
mass galaxies, the molecular gas is being consumed much 
more quickly than the atomic gas. 

(iv) The strongest dependence of id ep (H2) is on the spe- 
cific star formation rate of the galaxy, estimated using ei- 
ther the UV/optical light (corrected for extinction) or the 
4000 A break strength. This relation can be parametrized 
as tde P (H 2 ) = -0.44(log SFR/M* + 10.40) + 8.98, where 
tdep (H2) is in units of log yr _1 . There is also a strong de- 
pendence on stellar mass: tdep{i^2) = (0.36 ± 0.07)(logM* — 
10.70) + (9.03 ±0.99). 

(v) We found no significant correlation of the resiuals 
from these relations with any of the global galaxy properties 
we looked at. 

(vi) The tde P (H-2) versus sSFR relation extends smoothly 
with the same slope from the sequence of "normal" star- 
forming galaxies in our COLD GASS sample through to the 
population of more extreme starburst galaxies (LIRGs and 
ULIRGs), which have i dep (H 2 ) < 10 8 yr. 

(vii) Normal galaxies at z=l-2 are displaced with respect 
to the local galaxy population in the £d ep (H2) versus sSFR 



plane, in the sense that they have much longer molecular gas 
depletion times at a given value of sSFR. This is consistent 
with the fact that these galaxies are much more gas-rich. 



6 DISCUSSION 

In working towards the understanding of the star formation 
process, it probably makes sense to break the problem down 
into two parts: the conversion of atomic gas into molecular 
form, and the subsequent formation of stars from this cold 
dense phase. One way to address these issues is through 
resol ved studies of nea rby galaxy disks, on sub-kpc scales 
(e.g. iLerov et al.l l2008h or even down to t he scale of giant 
mole cular clouds (e.g. Gratier et al.l [201ol : ISchinnerer et al.l 
2010). Our approach is highly complementary. Only inte- 
grated measurements are available to us (total HI, H 2 and 
stellar masses, for example), but our COLD GASS sample 
is large, complete, unbiased, and covers a larger region of 
parameter space than previous studies. 

The central result of this paper is that the mean molec- 
ular gas depletion timescale (td ep (R2)) increases by a fac- 
tor of ~ 6 across the mass range sampled by COLD GASS 
galaxies, and that this variation correlates strongly with 
specific star formation rate, NUV— r colour and stellar 
mass. Previous surveys were unable to observe this vari- 
ation or correlation due to the limited dynamic range of 
their samples. It should be noted however that the non- 
universality of tde P (R2) is perhaps not very surprising. It 
is well known that gas depletion times are different in ex- 
treme starburst galaxies and merging systems in the local 
Univ erse, as well as submillimeter galaxies at high redshifts 
(e.g. iKennicuttl Il998bl: iRiechers et al.l [20071 : iBouche et"afl 
120071 : iBothwell et al.l 120101 ). In these galaxies, gas deple- 
tion times are significantly shorter, with the star formation 
surface density an order of magnitude larger at fixed gas 
surface density, com pared to the normal galaxy population 
jGenzel et alfoich . The novelty of our results lies in show- 
ing that tde P (R2) varies even within the population of nearby 
"normal" star forming galaxies, and in quantifying this vari- 
ation as a function of a range of fundamental parameters. 

Our work leaves open the intriguing question of what 
drives the variability in tde P (H2) across our sample. In the 
sma ller range o f galax ies e xplored in the TH INGS sam- 
ple, iBigiel et all (|2008l ') and ILerov et al.1 (120081 ) find a con- 
stant molecular gas depletion time. They do not observe 
any variation in this quantity with respect to galactocen- 
tric radius, Toomre Q, or any other galaxy propert y . This 
is consistent with the observation by iBolatto et al.1 (|2008h 
that molecular cloud properties appear essentially invari- 
ant across nearby galaxies with properties similar to those 
in the THINGS sample. In starburst galaxies, however, it 
is clear that tdep(H2) can be much small er even than the 
smallest values found in our sample (e.g. iGao fc Solomon! 
l2004h . and in such systems molecular cloud properties are 
measurably different than in normal galaxies (e.g. in M64, 
iRosolowskv fc Blitzl I2005T ). We therefore conjecture that 
the galaxies in our sample that show significantly reduced 
ide P (H2) are those undergoing minor starbursts, and that 
the weak correlation between td ep (H-2) and total stellar mass 
that we observe arises because smaller galaxies have more 
bursty star formation histories. In what follows, we investi- 
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gate this possible interpretation of the tdep (H2) relations, as 
well as two alternative explanations: quenching in high mass 
galaxies preventing the molecular gas from efficiently form- 
ing stars, and an increasing fraction of molecular gas not 
traced by the CO (1-0) observations at lower stellar masses. 

The first interpretation calls for enhanced star for- 
mation efficiency in low mass galaxies due to minor star- 
burst events. This picture is consistent with our result that 
tdep (H2) correlates inversely with specific star formation 
rate, and that the scatter in the relation is largest at low 
stellar masses. The starbursts in these galaxies would not 
be the large ones associated with major mergers, but would 
instead be weaker ones associated with distant tidal encoun- 
ters, variations in the intergalactic medium accretion rate, 
or secular processes within galactic disks. Any process capa- 
ble of driving gas toward the galactic center and raising the 
surface density there will suffice, since lower depletion times 
and altered molecular cloud properties are associated with 
galactic surface densities above ~ IO OMq pc~ 2 on both ob- 
serva tional and theoretical grounds l|Krumholz et alj|2006t 
120091) . 

The second interpretation for the variable td ep (H2) in- 
vokes the higher mass galaxies, and how star formation could 
be preferentiall y suppressed in thes e systems. Morphological 
quenching (e.g. Martig et alj I2009T ) and feedback scenarios 
which prevent the molecular gas from forming stars while at 
the same time not destroying the gas, would predict that a 
consequence of the quenching would be an increased reser- 
voir of molecular gas. The specific star formation rate would 
then inversely correlate with tdep(H2) as we observe, and we 
would expect to see some trend with galaxy morphology (e.g. 
concentration index, R90/R50) which is also seen in Figure 

As a third possible explanation, we ask whether the ob- 
served tdep (H2) trends may be explained by the fact that 
the CO (1-0) emission line may not be an accurate tracer 
of the molecular gas content of the low mass galaxies in 
our sample, which may have system atically lower metallic - 
ities than the h i gher m ass galaxies (|Tremonti et alj '2004). 
iKrumholz et al.l ()201 lh suggest that CO-inferred molecular 
gas measurements may begin to miss an increasing fraction 
of the total H2 at metallicities as high as half solar (see also 
iBolatto etal1l200Sl ; llerov et al.ll201lh . As discussed in Pa- 
per I, the COLD GASS sample is limited to galaxies with 
stellar masses greater than 1O 1O M0, so the global metallici- 
ties of our galaxies are all solar or greater, so this concern is 
not likely to apply. One might nevertheless still worry that 
metallicities in the outer regions of the galaxies are lower, 
so one might still not account for some fraction of the to- 
tal molecular content of the galaxy. If metallicity gradients 
were to correlate with for example stellar mass or stellar 
mass surface density, this might create the dependence of 
tdep (CO) on those quantities. However, our recent analysis 
of metallicity gradients using long-slit spectra obtained from 
MMT (Moran et al., in preparation) indicates that metal- 
licity gradients of the galaxies in our sample are mostly flat 
and do not correlate strongly with either quantity. 

The observed dependence of tdep (H2) on several param- 
eters such as stellar mass and specific star formation rate, is 
then likely caused by a combination of the first two effects, 
each acting at a different end of the COLD GASS stellar 
mass range. In high mass galaxies, star formation would 



preferentially be quenched by internal processes, while in 
low mass galaxies, star formation would preferentially be 
enhanced through mild starburst events. 

The final mystery in this picture is why tdep (HI) should 
vary less than tdep (H2). We note that the invariance of 
^dep (HI) must break down in sufficiently strong starbursts 
- for exa mple Arp 220 has a star formation rate of ~ 50 
Mq yr _1 (Dowries & Solomon 1998), but an HI mass below 



4 x 10 9 M Q (|Baan et al.lll987l ). giving t dep (HI)< 100 Myr. 



However, galaxies like Arp 220 are rare, and we almost cer- 
tainly do not have any galaxies quite that extreme in our 
sample. For the weaker starbursts we do include, it seems 
plausible that the majority of the HI resides in a quiescent 
outer disk that is largely unaffected by whatever process 
produces the starburst. It indeed appears that the conver- 
sion of atomic to molecular gas is, in some situations, a bot- 
tleneck in the star formation process. For example, at low 
stellar mass surface densities (log^i, < 8.7, see Figure [3p), 
all galaxies have td ep (H2) shorter than the sample mean of 1 
Gyr. These galaxies are processing their molecular gas into 
stars faster than they are consuming their atomic gas, due 
to some internal (e.g. bar instabilities) or external (e.g. tidal 
interactions) processes triggering minor starburst events, as 
we speculated above. In these galaxies, HI is likely to be ac- 
cumulating in regions where the conversion between atomic 
and molecular gas is less efficient and where this gas cannot 
participate in the the star formation process, for example in 
the outer disks of galaxies. 

Why should HI accumulate more effectively in low mass 
galaxies than in high mass galaxies? One possibility is that 
gas from the external environment is able to accrete more 
efficiently in such systems. At a halo mass of greater than 
few xlO n M0, gas is expected to shock-heat to the virial 
temp erature of the ha lo and cooling times become longer 
(e.g. iKeres et al.ll2005h . Once the gas has been accreted, 
the production of molecular gas may rely of gas flows that 
bring gas from the outer disks to the central regions of the 
galaxy, where it is able to reach high enough densities to 
be shielded from the ambient UV radiation field and form 
molecules. If such gas flo ws are dynamically driven (e.g. 
IChakrabarti fc Blitj Boblh . one might expect them to oc- 
cur more frequently in high mass galaxies, which are located 
in more massive halos and more crowded environments. In 
addition, feedback processes from AGN may disrupt gas ac- 
cretion in the most massive dark matter halos. 

It is also interesting to speculate on the issue of gas ac- 
cumulation in high redshift galaxies. We have seen that the 
molecular gas depletion times in normal galaxies at z — 1 — 2 
are nearly an order of magnitude longer than in local galax- 
ies with the same specific star formation rates. In nearby 
starburst galaxies, the very high values of sSFR are achieved 
because dynamical disturbances act to compress the avail- 
able interstellar medium into a very small volume, in which 
the filling factor of GMCs and young stars is presumably 
very high. Although we have not yet been able to demon- 
strate this conclusively, our current best guess is that the 
relation we uncovered between tdep (H2) and sSFR is sim- 
ply a reflection of increased dynamical stirring of one sort 
or another, as one proceeds to more strongly star-forming 
galaxies. At the extreme high sSFR end of the sequence, 
clear evidence for recent merging is always seen, but at the 
low sSFR end, the diagnostics of the dynamical effects that 
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are at work may be much more subtle. We are currently in 
the process of analyzing whether bar-driven inflows, for ex- 
ample, might account for some of the trend internal to our 
COLD GASS sample. 

In high redshift galaxies, however, the high values of 
sSFR are achieved simply because there is a lot of molecular 
gas. This leads us to speculate that at high redshifts the 
star formation bottleneck is no longer in the conversion of 
atomic to molecular gas. Indeed, the flat dependence of Qhi 
with redshift inferred from the abundance of damped Lyman 
alpha systems in high redshift quasar spectra would appear 
to rule out th e presence of large HI reservoirs in high redshift 
galaxies (e.g. iHopkins et all [2005). Although gas accretion 
rates are expected to be much larger at high redshifts, most 
of the gas must be transformed to molecular form almost 
immediately. The star formation bottleneck in high redshift 
galaxies may then have more to do with the formation of 
the proto-stars themselves. 
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APPENDIX A: STAR FORMATION RATES 

We use a spectral energy distribution (SED) fitting tech- 
nique to derive the star formation rate (SFR) for each galaxy 
in the COLD CASS sample. The lBruzual fc Charlotl (f2003T ) 
population synthesis code is used to generate model galax- 
ies with a range of ages, star formation histories, metallici- 
ties and dust attenuations and to create a library of model 
SEDs in 7 bands (FUV, NUV, u, g, r, i and z). For each 
galaxy, we evaluate the goodness-of-fit of each model SED 
and then estimate its SFR and dust attenuation factor as 
the x 2 -weighted average of the model p arameters. 

This technique is very similar to ISalim et alj (|2007T ) 
in the way the metallicity, age and star formation history 
(SFH) priors are set. There are, however, two main differ- 
ences compared to their methodology: (1) we only use an 
exponentialy declining SFH without bursts, because mas- 
sive galaxies such as those in the COLD GASS sample do 
not tend to have bursty SFHs, and (2) we use an extinction 
coefficient (Av) that depends on NUV— r colour, after cal- 
ibration using a reference sample with direct measurement 
of the SFR through UV and infrared measurements. The 
technique and the calibration of Av are explained in more 
detail below. 

As pointed out by ISalim et al.1 (120071 ) and IWang et all 
(2U l it), the derived dust attenuation (and thus SFR) is sen- 
sitive to the assumed prior distribution of Av- A direct es- 
timate of Av in a galaxy is possible only if we know its UV 
through far-IR SED. This is not the case for our sample, 
so we have assumed a different prior Av distribution for 
each individual galaxy, constrained by an indirect estimate 
of Av in the galaxy. Two different methods of estimating the 
Av prior are considered: the measured Balmer decrement 
from the SDSS s pectrum, and a com bination of D4000 and 
NUV-r colour (| Johnson et al.l 120071 ). The former method 
measure s the extinction exp erienced by stars located in HII 
regions ( Calze ttTet al.|[l994l ) in the central 3" of the galaxy, 
while the latter is a rough estimate of the attenuation of the 
ligh t from stars with ages of a few hundred million years (see 
also I Schiminovich et al-lfeoiol ). We denote the Av measured 
from these two methods as Ay(B) and Ay(N), respectively. 
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Figure Al. Difference in the SFRs measured using the SED-fitting technique and a combination of UV and IR photometry 
(ASFR = log(SFR SBD ) - log(SFR t0 i)). The values of ASFR are plotted as a function of NUV-r colour (top panel), and the 
colour distribution of the COLD GASS sample is shown in the bottom panel (in gray the complete sample, in blue the galaxies 
with CO detection). Over the colour range where CO detections are possible (NUV— r< 5), the scatter between the two SFR 
measurements is only 0.24 dex, with no systematic offset. 



We make use of the GMACS (Gala xy Multiwave- 
lenth Atlas From Combined Surveys) sample ijjohnson et al.l 
120071 ) to tune our choice of the assumed prior distribution 
of Av in the model library, and to test the robustness of the 
SFR derived from the SED fitting. The 70^tm fluxes are con- 
verted into total infrared luminosities (Ltir), taking into 
account a correction facto r from the ratio of t he observed 
fluxes at 24/xm and 8um dBoquien et al.ll2010l ). We apply 
the dust m odel from Meurer et al.l fl_999|) and the extinction 
curve from lCalzetti et ah 1 1| 19941) and calculate Ay from the 
ratio of Lt ir over the lumino sity in the FUV (j4y(IRX)). 
We use the iKennicuttl |l998al ) relation to convert Ltir to 
SFR (SFRt/jj) and divide it by a factor of 1.7 to account for 
the difference between Salpeter IMF and Chabrier IMF. We 
calculate SFR from o bscured FUV lumin osity (SFRfc/v) us- 
ing the formula from lSalim et all (|2007l ). The total star for- 
mation rate for the GMACS galaxies is therefore obtained 
as SFR to t=SFRT/ii+SFR(7y , with the latter not corrected 
for dust extinction. 



Different prescriptions for Av based on the SDSS data 
are tested against the directly measured values of Ay{lRX). 
We find that different techniques to measure for Ay min- 
imize best the scatter with ylv'(IRX) in different NUV— r 
colour intervals, and therefore adopt the following prescrip- 
tion: 



Ay (best) 



A V (B) - 0.052(NUV - r) 
0.5{Av(B) + A V (N)) 
A V (N) 
1.5 



if NUV-r< 2.5 
if 2.5 <NUV-r<3.5 
if 3.5 <NUV-r< 5.0 
if NUV-r> 5. 



(Al) 



Therefore, when we estimate the SFR for an individual 
galaxy, we adopt a Gaussian prior distrubtion of Av based 
on its observed colour that peaks at Av{best) as given in 
Equation IA1 1 and has a width of a = 0.15. We compare SFR 
from the SED fitting method with the "true" SFRtot mea- 
sured from the UV and IR photometry in Figure IA1I The 
two sets of SFRs are very consistent when NUV— r< 5, with 
a scatter of 0.24 dex in their log difference. The scatter over 
the full colour distribution is 0.37 dex, because of the low 
S/N in the GALEX images for the redder galaxies. This, 
however, has no influence on the results of this paper, since 
all COLD GASS galaxies detected in CO have NUV-r< 5 
where the SED-fitting technique produces reliable SFR mea- 
surements. 

Another approach to measure SFRs for the COLD 
GASS galaxies would have been to correct the GALEX FUX 
flux using ^^(best) to infer the total SFR (SFRFc/v,cor), in- 
stead of relying on SED fitti ng using a mo d el lib rary. This 
was the approach taken in I Johnson et al.l (|2007l ). and we 
test it on the GMACS sample. The scatter in the log dif- 
ference of SFRtot and SFKfuv,cot is 0.29 dex for galaxies 
with NUV— r< 5. Over the same colour range, the scatter 
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between SFR to t and SFRsbd is 0.24, indicating that infor- 
mation is gained by using the full 7 optical and UV bands 
from SDSS and GALEX in the determination of the SFRs, 
if the dust attenuation prior distribution can be well con- 
strained. 



APPENDIX B: RESIDUALS OF THE 
MOLECULAR GAS DEPLETION TIME 
RELATIONS 

As explained in Section 13.11 and shown in Figure [3] the 
molecular gas depletion time (tdep(H2)) correlates best with 
specific star formation rate (sSFR). To investigate what 
other physical parameters may be at play in setting this re- 
lation, we examine here the residuals of the sSFR-t<j e j,(H2) 
relation as a function of the other global parameters con- 
sidered in this study. The (data - fit prediction) residuals 
are shown in Figure [BTI where a positive residual therefore 
indicates a longer depletion time than average at fixed sSFR. 

The only two parameters upon which the residuals de- 
pend are the measures of smoothness and asymmetry (Fig- 
ure lBll h and i). However, the trends are driven mostly by a 
small number of points at both low and high values of Smth 
and Asym, and the amplitude of the trend seen in the binned 
values is smaller than the overall calibration accuracy of the 
star formation rates (Appendix |X| and the uncertainty on 
the molecular gas masses. Furthermore, if we were instead 
to show the residuals of the stellar mass-^^H^) relation, 
the trend seen with Smth and Asym is even weaker. We 
therefore conclude that there is no significant third parame- 
ter dependence among the quantities we have examined. We 
will exmine this issue in more detail in a later paper. 
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Figure Bl. Residuals of the sSFR-i^ej, (Hz) relation of Figure [3};, as a function of 12 different global measurements for the 
COLD GASS galaxies with secure (filled circles) and tentative (open circles) CO detections. The red symbols show the mean 
value of the residuals in equally-populated bins of the a;— axis parameter. The correlation coefficient (r) of the data in each 
plane as well as the scatter around the best fit relation (tr) are given in each case. 
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